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ABSTRACT 

We report on the SDO/AIA and Hinode/ElS observations of a transient coronal loop. The loop 
brightens up in the same location after the disappearance of an arcade formed during a B8.9-class 
microflare three hours earlier. EIS captures this loop during its brightening phase as observed in most 
of the AIA filters. We use the AIA data to study the evolution of the loop, as well as to perform 
the DEM diagnostics as a function of ft. Fe XI-Fe XIII lines observed by EIS are used to perform 
the diagnostics of electron density and subsequently the diagnostics of ft. Using ratios involving the 
Fe XI 257.772A selfblend, we diagnose ft <2, i.e., an extremely non-Maxwellian distribution. Using 
the predicted Fe line intensities derived from the DEMs as a function of ft, we show that, with 
decreasing ft, all combinations of ratios of line intensities converge to the observed values, confirming 
the diagnosed ft <2. These results represent the first positive diagnostics of ft-distributions in the 
solar corona despite the limitations imposed by calibration uncertainties. 

Keywords: Techniques: spectroscopy - Radiation mechanisms: non-thermal - Sun: corona - Sun: UV 
radiation - Sun: X-rays, gamma rays 
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The distance and nature of the stellar coronae with 
temperatures of up to several million Kelvin together 
with the proximity of the stars themselves make direct 
probing of these environments impossible at present. In 
the absence of in-situ measurements, the emitted spec¬ 
trum remains the only source of information about the 
physical conditions in these emitting media. Most of 
the coronal radiation is emitted in the X-ray, extreme 
ultraviolet (EUV) and ultraviolet parts of the electro- 
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magnetic spectrum, requiring space-borne observatories. 
This poses further problems with limi ted availability of 
the data, calibration a nd its stability (Be nMoussa et a U 
l2Q13t iDel Zannall2Q13al ) , as well as the unavoidable trade¬ 
offs between the spatial, temporal, and spectral resolu¬ 
tions and coverages. 

In many instances, it is advantageous to combine imag¬ 
ing and spectroscopic observations (e.g.. iSchmelz et al.l 
120091: I Warren et al.l 120121: IDel Zann4 l2013bh . Imaging 
observations employing narrow-band EUV filters, such as 
those made bv the Atmospheric Imagin g Assembly (AIA, 
iLemen et al.l 1201 3 iBoerner et al. 1201211 onboard the So¬ 


lar Dynamics Observatory fSEO lPesnell et al.ll2012lh of¬ 
fer high spatial and temporal resolution. These are com¬ 
plemented with spectroscopic observations made in sim¬ 
ilar wavelength ranges, such as th ose performed by the 
EUV Imaging Spectrograph ( EIS. iCulhane et~ah 120071) 
onboard the Hinode spacecraft (|Kosugi et al.ll2007 ). Tra¬ 
ditionally, the vast majority of both the imaging and 
spectroscopic data are analyzed and modelled under the 
assumption of the local equilibrium, i.e., Maxwellian dis¬ 
tribution of particle energies. This is because calculation 
of the synthetic spectra requires integration of many in¬ 
dividual ionization, recombination, excitation and deex¬ 
citation cross-sections over the (unknown) distribution 
function in order to obtain rates of these processes and 
finally the corresponding emissivities at individual wave¬ 
lengths. The assumption of Maxwellian distribution af¬ 
fords easy calculation of the synthetic spectra , e.g., usinfi 
the CHIANTI atomic d atabase and software (|Dere et al 

ITflflTt ILandi ^~7Ill20T3h . 

However, this assumption is incorrect if there are 
correlations between the particles in the system. 
Such correlations can be ind uced by any lo ng-range 
interactions in the system (iCollierl 12004 iLeubnerl 
12004 iLivadiotis fc McComasI I2009L l2010l I2013D . e.g., 
parti cle acceleration due to magnetic reconnection 
(e.g., Zharkova et al. 2011; Petkaki & MacKin nonl 2011 

~ inr 


Stanier et al. 2012; Cargill et al. 2012; Burge et al. 2012 
2014 iGordovskvv et al.l I2Q13L 120141). s hocks , or wave- 
particle interactions (e.g., IVocks et al.l 120081 ). In such 
cases, the particle distribution will depart from the 
Maxwellian one, and will likely exhibit an enhanced 
high-energy tail. Furthermore, turbulence with the dif¬ 
fusion coefficient inversely proportional to particle ve¬ 
locity will also lead to the appearance of the non- 
Maxwellia n distributions with characteristic high-energy 
tails fe.g.. lHasegawa et al.lIl985t lLaming fc Leoril 120071: 
iBian et al J I20l4h We note that the collision cross¬ 
sect ion^caJ^s_wrthE -2 , where E is the particle en¬ 
ergy ([Mever-Vernetl I2007D . This leads to the be¬ 
haviour of the collision frequency as E -3 / 2 , i.e., the 
high-energy tail is difficult to equilibrate. Therefore, 
the currently fav oured theories of nanoflare heating of 
the solar corona (iKhmchukl 120061: iKlimchuk et al.l [20101: 


Tripathi et al. 2010; Bradshaw et al. 2012; Winebarger 
20121: Viall & Klimchuk 2011a.b, 2013) should afford sit- 


uations for departures from th e Maxwellian distribution. 
iScudder fc Karimabadil ( 201 3) argue that the particle 
distribution in the solar and stellar coronae above l.O5i?0 
is strongly non-Maxwellian. 

Observational clues that the solar corona could 
be non-Maxwellian come from the in-situ detection 


of the non-Maxwe l lian /^- distributions in t he sol ar 
wind (ICollier et al.l 119961: iMaksimovic et a hi Il997allbl : 
IZougane lis 2008; Le Chat et al.ll20lTl) . Spectroscopic ev¬ 
idence for the presence of ^-distributions in the transi- 
tion region was found from a ctive region Si III spectra 
(|Dzifcdkova fc Kulinoval 1201 lh . The /^-distributions are 
characterized by a high-energy power-law tail (Sect. [2]), 
with the power-law index being given by the value of 
— (ft-f 1/2). In the solar wind, the typically detected 
values are ft >2.5, while the detection in the transition 
region yielded ft ~ 7 from the active region Si III spec¬ 
trum. Presence of high-energy electrons in the transition 
region at the base of coronal l o ops h ave recently been 
established also by iTesta et al.l (12014 by analyzing the 
transition -region Si IV emission ob served by the IRIS in¬ 
strument (|De Pontieu et al.ll2014D . 

Despite these detections in the solar wind and the 
transition region, no direct, unambiguous detection of 
the non-Maxwellian distribution in the coronal spectra 
have been mad e to date. At t empts at doing so were 
made e. g. by IFeldman et al.l (120071). iRalchenko et al. 
2007D . lHannah et al.l j 201 Oil . iDzifcakova fc Kulinova 


(2010), and IMackoviak et all (120131 ). IFeldman et al 


( 2007 ) assumed a bi-Maxwellian distribution, with the 
temperature of the second Maxwellian chosen to be 

10 MK, and argued that no such second Maxwellian 
is necessary to explain the observed spectra of He-like 
lines. This analysis was however limited and did not 
include the effe c ts of ft-distributions on the spectra. 
IRalchenko et al.l (2007) showed that the quiet-Sun Si, 
Ca, and Ar spectra are consistent with a bi-Maxwellian 
distribution, where the second Maxwellian contains only 
a small fraction of particles, of the order of several per 
cent. Its temperatures range from 2.3 MK (300 eV) to 
7.7MK (lOOOeV), with the fraction of particles being 
at most 5-7%, and 1% for these temper atures, respec¬ 
tively. No diagnostics of ft was performed. lHannah et al.l 

11 used off-limb qui et-Sun observati ons performed by 
the RHESSI spacecraft (|Lin et al.lll2002 ) to constrain the 
emission measures corresponding to individual power- 
law and ft-distributions. However, the constraints ob¬ 
tained on the emissi on measures were rather large even 
for small values o f ft. [D zifcak ova fc Kulinoval j2010) and 
IMackoviak et al ] (I 20 T 1 studied the possibilities of di¬ 
agnosing ft using lines observed by EIS. Despite having 
found several combinations of line ratios sensitive to ft, 
mostly involving neighbouring ionization stages, unam¬ 
biguous detection were difficult due to atomic data un¬ 
certainties and possible multi-thermal effects along the 
line-of-sight that were n ot accounted for. Subsequently, 
IMackoviak et al.l ( 2014 ) showed that the techniques to 
obtain the differential emission measures work also for 
the ft-distributions, and studied the influence of ft on 
such analyses. 

In this paper, we use the imaging and spectroscopic 
observations performed by SDO/AIA and Hinode/EIS 
in conjunction with the latest instrument calibration, 
atomic datasets, and differential emission measure tech¬ 
niques to analyze a transient loop observed within an 
active region core. Spectral synthesis for the non- 
Maxwellian ft-distribution is briefly discussed in Sect. O 
Analysis of the AIA imaging observations is performed 
in Sect. [3] In Sect. |4j ratios of spectral lines observed 
by EIS are analyzed to obtain the electron density and 
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log (T/K) = 6.00 



Figure 1. The /^-distributions with k = 2 : 3, 5, 10, 25 and the 
Maxwellian distribution plotted for log(T/K) = 6.0. Colors and 
linestyles denote the different values of k. 

A color version of this image is available in the online journal. 


Fe ionization equilibrium 



Figure 2. Ionization equilibrium for Fe X I-Fe XIII for the /re¬ 
distrib utions, according to the calculations of Dzifcakova &; Dudfkl 
(f20T3h . 

A color version of this image is available in the online journal. 

ft in the loop. Influence of the atomic data uncertainties 
on the diagnostics are discussed in Sect. [5] The results 
are summarized and discussed in Sect. [6j 

2. SYNTHETIC SPECTRA FOR THE /^-DISTRIBUTIONS 
2.1. The Non-Maxwellian k- distributions 

The /^-distribution of electron energies (Fig. [T|) is de¬ 
fined as a two-parametric distributi on with parameters 
T G ( 0, +oo) an d ft £ (3/2, +co) (e.g ., lOwocki fe Scndderl 
119831 : ILivadiotis fc McComasI 120091) 

9 F’^/^dF 

UE)dE = - - — +r, (1) 

V 1 ^ («—3/2)/c b T J 

where the A K = T(ft + 1)/(T(ft — 1/2 )(k — 3/2) 3 / 2 ) is the 
normalization constant and &b = 1-38 xlO -16 ergs -1 is 
the Boltzmann constant. The Maxwellian distribution at 
a given T is recovered for ft—>> 00 . Maximum departure 
from the Maxwellian occurs for ft —3/2. 

The mean energy (E) = 3&bT/ 2 of a ft- 

distribution does not depend on ft. Because 


of this, T can be defined as the temperature. 
iMever-Vernet et"ahl (1199511 . ILivadiotis fe McComasI 

(12009ft and ILivadiotis fe McComasI (|2010l ) showed 
that T indeed has the same physical meaning for the 
ft-distributions as the (kinetic) temperature for the 
Maxwellian distribution, and that it also corresponds 
to the definition of physical temperature in the frame¬ 
work of the gener alized Tsallis statistical mechanics 
(jTsallisI fl988l 2009 ). This permitted these authors to 
generalize of the zero-th law of thermodynamics for the 
ft-distributions, and it also permits, e.g., the definition of 
electron kinetic pressure p = n e k^T in the usual manner. 

It is straightforward to see from Eq. O that in the 
high-energy limit, the ft-distribution approaches a power- 
law with the index of —(ft + 1/2). On the other hand, 
ILivadiotis fe McComasI (j2009l ) showed that, in the low- 
energy limit, the ft-distribution behaves as a Maxwellian 
with Tm = T (ft — 3/ 2) / (ft + 1), and scaled to an appro¬ 
priate constant (see [Dzifcako va et all l2015h . Therefore, 
the ft-distribution can be thought of as a Maxwellian core 
(at a lower temperature) with a power-law tail. 


2.2. Line Intensities for the ft- Distributions 

The emissivity Sji of an optically thin spectral line 
arising due to a transition j -A i, j > i, in a &- 
times ionized ion of the element X is given by (e .g., 
iMason fe Monsignori Fossilll994 iPhillips et alJl2QQ8h 


£ 


ji ~ 


— A r,(X+ k ' \ — Li n A +fc ) n(X+ fc ) 

A ji " 1 , ) He n{x+k) n{x) 

AxGx,ji(T , n e , ft)n e nH , 


A x n e nn 


( 2 ) 


where h~6.62 x 10 -27 ergs is the Planck constant, 
c~ 3 xlO 10 cms -1 represents the speed of light, Xji is 
the wavelength of the resulting spectral line, Aji the cor¬ 
responding Einstein coefficient for spontaneous emission, 
n(X+ k ) is the density of the ion +k with electron in the 


excited upper level j, n(X +k ) the total density of ion 
+fc, n(X) = n x the total density of element X whose 
abundance is A x , ur the hydrogen density, and n e the 
electron density. The function G X ji(T , n e , ft) is the con¬ 
tribution function for the line Xji. The G X ji(T, n e , ft) 


is a function of ft due to the dependence of the in¬ 
dividual collisional processes of ionization, recombina¬ 
tion, excitation and deexcitation on ft (e.g.. iDzifcakoval 


19^1 lDzifr,akav^l200l 2006a 

Wannawichian et al.l 

200 IJ 

Dzifcakova & Masonl 20081: 

Dzifcakova fe Dudfk 

2013; 

Dudfk et al.ll2014a.b: Dzifcakova et al.ll2015D. 


The intensity Iji of the spectral line is then given by 
the integral of emissivity along a line of sight l 


or 


Iji = J A x G X ji(T,n e , K)n e nudl , 

Iji = J A x G X ji(T, 7i e , ft)DEM^(T)dT, 


( 3 ) 

( 4 ) 


where the quantity EM = J n e nudl is the plasma emis¬ 
sion measure, and DEM^(T) = n e 7iRd//dT is the dif¬ 
ferential emissi on measure, generaliz ed for the ft- 
distributions bv iMackoviak et all (j2014l h 
The line intensities for ft-distributions are evalu¬ 
ated using the ionization equilibrium calculations for 
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/^-distributions (|Dzifcakova fc Dudikl I2013D . Relative 
level popul ation is obtai ned us in g the approximative 
metho d o f iDzifcakoval (1200614). IDzifcakova fc Masonl 
(j2008l ) and IDzifcakova et al.l (j2Q15h based on atomic data 
corresponding to the CHIANTI da tabase, version 7.1 
(|Dere et al.l 119971: lLandi et al.l 120131) . The accuracy of 
this method for allowed tr ansitions is less than 10 % 
(|Dzifcakova & Masonl I2008D . For the Fe IX-Fe XIII, 
which are the ions used for diagnostics in Sect. 0J 
the line intensities are ob tained by di rect i ntegration of 
the collision strengths (|Dudfk et al.ll2014bh . The colli¬ 
sion strengths fo r these calc u lation s are state-of - the-ar t 

” (Ini, 
( 2om 
B . 


Del Zanna et al. 


12012 bl 
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IDel Zanna & Storev 

Del Zanna et al. 
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3. SDO/AIA OBSERVATIONS OF A TRANSIENT CORONAL 
LOOP: IMAGING 

The Atmosphe ri c Imaging Assem bly (AIA, 
iLemen et al.l 120121: iBoerner et al.l I2012D on board 
NASA’s S o lar Dynamics Observatory (SDO, 
iPesnell et al.l I2012D mission consists of 4 identical, 
normal-incidence two-channel telescopes with a diame¬ 
ter of 20 cm. AIA provides multiple, near-simultaneous 
full-Sun images with high temporal (12 s) and spatial res¬ 
olution (1.5", pixel size 0.6"). AIA images of the Sun are 
taken in 10 filters, 7 of which are centered on EUV wave- 
lengths (94A, 13lA, 17lA, 193A, 21lA, 304A,and 335A), 
and 3 on UV or visible wavelengths (1600A, 1700A, 
and 4500A). The EUV filters are centered on some 
of the strongest Fe lines in the solar EUV spectrum. 
However, other emission lines originating at different 
temperatures are present within each filter bandpass. 
This makes the temperature respo nses of the AIA EUV 
filters highly multi-thermal (e.g., IQ ’Dwyer et al.l 120101: 
iDel Zanna et al.ll2011bl: iSchmelz et al.ll2013l) . 

Compared to the Maxwellian distribution, the AIA re¬ 
sponses for /^-distributi ons are even more multi-thermal 
(|Dzifcakova et al.lI2015L Fig. 8 therein). This is mostly 
caused by the behav iour of the ionization equilibrium 
(|Dzifcakova fc Dudfkll2013lh with the relative ion abun¬ 
dances of individual ions becoming wider and flatter with 
decreasing n (see also Fig. [ 2 ]). Peaks of the individ¬ 
ual filter responses can also be shifted towards higher 
T, especially for low k = 2-3 (jDzifcakova et all 120121: 
IDzifcakova et al.ll2015l) . Typically, in-depth understand¬ 
ing of the AIA observations requires d ifferential emission 
measu re ( DEM) analysis, see , e.g., lHannah fc Kontarl 
(| 2012 l) and lHannah fc Kontarl (|2013l) . 

In the following, we report on the AIA observations of 
a transient coronal loop and its spatial and temporal re¬ 
lation to the solar microflaring activity. DEM analysis of 
the AIA data is performed in Sects. 13.31 and its influence 
on spectroscopic diagnostics is studied in Sect. 14.31 


3.1. Relation of the Transient Loop to Microflaring 
Activity 

The transient loop studied here was observed in the 
active region (hereafter, AR) NO A A 11704. This active 
region was a small bipolar AR, of Hale class /3/a, char¬ 
acterized by two leading negative sunspots and dispersed 
plages of both polarities. On 2013 March 30, a compact 


B 8 .9-class microflare was observed within this AR, peak¬ 
ing at about 10:24UT (Fig. [3j top left). A weaker B4.8 
microflare was observed to peak at 13:21 UT in a differ¬ 
ent AR NOAA 11708. Perhaps by coincidence, during 
the same time, we observe a transient loop within AR 
11704 in the same place as the previous B8.9 microflare 
(Fig. 12]). 

Figure |3] and Movie 1 (online) show the location of the 
B8.9 microflare and the transient loop within the AR 
11704 in AIA 171A channel. The 171A channel is chosen 
since it provides the best representation of the morphol¬ 
ogy of all features. The B8.9 flare was related to a failed 
eruption of a long intermediate filament. One end of the 
filament lies within AR 11704, with the filament extend¬ 
ing further away in the North-East direction (Fig. [3]). 
The activation of the filament starts at about 10:00 as a 
series of brightenings within the filament and its imme¬ 
diate neighborhood, subsequently leading to apparently 
untwisting motions of bright threads. The filament does 
not erupt however, perhaps due to the overlying field 
seen e.g. in the AIA 193A band. 

The failed eruption is accompanied appearance of a 
flare arcade within AR 11704. The arcade first appears in 
the AIA 13lA band (Fig. [3j top right) and subsequently 
in cooler AIA channels. The arcade is still visible in 
17lA an hour later at 11:24 UT and subsequently fades. 
The fading of the arcade is accompanied by the appear¬ 
ance of falling blobs along the arcade in 304A (Movie 
1) that are visible until about 12:30 UT. After this time, 
evolving loops are still observed in the same location. 
The evolution of these loops (Movie 1) resembles slip - 
ping reconnection (e.g JAulanier et al.ll2006[f2007l 120121) . 
although this is ambiguous given the presence of many 
individual loops and other emitting structures. The loop 
emission at the loci of the flare arcade completely fades 
out of the AIA 17lA channel at about 12:53 UT. How¬ 
ever, at approximately 13:01 UT, a transie nt c oronal loop 
starts brightening up once more (Sect. 13.2)) . preceded 
again by activity in the filament (Fig.[3j bottom left). No 
restructuring of the large-scale magnetic configuration of 
the AR 11704 is noticeable during the period studied here 
(09:00-15:00UT, Fig. 0. 

3.2. Multi-wavelength Evolution of the Loop 

We now focus on the evolution of the transient loop. 
The multi-wavelength evolution of the EUV emission is 
shown in Figure 0] with a cadence of 20 minutes. The 
emission of the transient loop is observed in AIA 171A for 
nearly 2 hours (see also online Movie 1). The images 
in Fig. 0] are chosen to include the time of 13:19 UT, 
during which the EIS spectrometer observes a part of 
the transient loop (Sect. 0j). 

Upon careful examination, we find spatial misalign¬ 
ments of unknown origin between various AIA bands, 
especially in Solar Y direction. These misalignments 
are corrected for manually by matching the positions 
of the low-lying moss emission. In this correction, the 
AIA 171A and 193A bands are taken as the reference 
ones, since no detecable misalignment between these two 
bands is found. The AIA 211 A is corrected by shifting it 
by AY = — 1 . 6 px (lpx = 0.6") to match the 193A band, 
which has similar emission morphology. Subsequently, 
the 335A band is corrected by shifting it by AY = 2.13 px 
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Figure 3. Evolution of the two microflares observed by AIA 171A. Top left : GOES X-ray flux at 1-8A. Red-dashed vertical lines denote 
the times of the two microflares, full brown lines denote the times of AIA 171A snapshots, and the orange-dashed lines denote the times of 
snapshots in Fig. [4] Inset in top right : AIA 131A 1-minute average image showing the arcade of during the B8.9 microflare. The location 
of the inset corresponds to the box shown. 


A color version of this image is available in the online journal. 

to match the corrected 211A band, again due to sim¬ 
ilar emission morphology. The 13lA band is shifted 
AY = 1.5px to match the 17lA band. No correction 
is necessary for the 94A band. After these corrections, 
the location of the moss emission is the same through¬ 
out the AIA bands, and the spatial correspondence of 
loops observed in different bands is substantially im¬ 


proved as well. We also note that the AIA data shown 
in Fig. 2] are corrected for solar differential rotation using 
the drot-map routine available within Solar Soft under 

IDL. 

Before the loop appears in AIA 171 A, it is visible 
in AIA 335A, 21lA, and 193A (Fig. IH top row). The 
transient loop belongs to the series of evolving loops al- 
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Figure 4. Multi-wavelength imaging observations of the loop evolution made by SDO/ AIA. The time cadence shown is 20 min. AIA 
data are averaged over 1-minute intervals to increase the signal-to-noise especially in the 335A and 13lA bands. The images are scaled 
logarithmically. The field of view shown corresponds to the box in Fig. [3] Small black and green boxes in the second row denote the 
selected portion of the loop and background, respectively. The white box shows the field of view of a portion of the EIS raster shown in 
Fig. [8] 

A color version of this image is available in the online journal. 
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Figure 5. Top : Background-subtracted intensity lightcurves for 
the small black box shown in Fig.0] The lightcurves are normalized 
to maximum intensity in the time window shown. The maximum 
intensities are listed for each filter. Individual colors stand for 
individual AIA filters. Bottom: Time-distance plots in AIA 131A, 
17lA, and 193A along a slit placed through the black and green box 
at Solar Y = 392.7 // . Black and green horizontal lines correspond 
to the Solar X extensions of the respective boxes shown in Fig. 0] 
where black stands for the loop and green for the background. The 
vertical dashed black line corresponds to 13:19 UT (Sect. |4j). 

A color version of this image is available in the online journal. 


ready present at the same location (Movie 1). At about 
13:01 UT (Fig.®, the loop is discernible as a brighten¬ 
ing structure in all AIA EUV channels. In the vicinity 
of its right footpoint at X=235”, Y=370”, we observe 
heightened activity (see Movie 1), with the appearance 
of many bright, short, short-lived closed loops and jet¬ 
ting activity, likely indicating ongoing magnetic recon¬ 
nection. The transient loop continues to brighten, and 
at 13:19 UT, corresponding to the time of the EIS obser¬ 
vation, the loop is already a prominent emitting structure 
(Fig. 01 second row). The temporal evolution of the loop 
emission is shown in Fig. [5j The top panel of this fig¬ 
ure shows background-subtracted lightcurves in a small 


black box of 4x5 AIA pixels at the position of the loop. 
The background is represented by the green box of 4x5 
AIA pixels located to the right of the black box (Fig. 01 
second row). Its location is chosen so that its intensity 
is not contaminated by the jet emission (Fig. [5]). 

The loop emission first peaks in the 335A filter at about 
13:20 UT, followed by 21lA and 193A at 13:26 UT, and 
subsequently by 94A at 13:27 UT, 171A at 13:29 UT, and 
13lA at 13:32 UT that order. I.e., the emission first 
peaks in Fe XVI and subsequently cools down to Fe VIII 
observed in 13lA. This is important, since these results 
show that no hot flare-like emission is present. If it were, 
it would be detected as Fe XXI and F e XVIII emission 
in the AIA 131A and 94A bands (e.g., iDel Zanna et al.l 
i2ona iPetkaki et al.ll2012h prior to the maximum of the 
335A lightcurve. The 94A emission is observed to be very 
weak, of the order of several DN s -1 px -1 , suggesting ab¬ 
sence of strong, dense plasma emitting in Fe XVIII. The 
131A emission closely follows the 171 A, showing that it 
is dominated by Fe VIII. This conclusion is confirmed 
using the time-distance plots constructed as a function 
of Solar X using a slit placed at Solar Y = 392.7", i.e., 
through the centres of the black and green boxes (Fig. 01 
second row). 

We note that the lightcurves show multiple peaks. This 
is most likely due to the multi-thermal nature of the 
AIA bands. In principle, DEM analysis can help identify 
the contributions to individual peaks, if the distribution 
function (i.e., value of k) is known. Since n cannot be 
determined from AIA data alone (Sect. 13.3 1) . we do not 
perform this analysis. 

We also note that the loop evolution is accompanied by 
presence of multiple jets originating near the loop’s right 
(western) footpoint (Movie 1). With increasing height, 
these jets diverge from the loop and thus do not contami¬ 
nate the loop emission. In Fig.[5j bottom , the jets appear 
as a series of fast-moving, short-lived brightenings near 
the position of the loop (black horizontal lines). The 
jetting activity starts at about 13:14 UT and extends be¬ 
yond the lifetime of the transient loop. 

With increasing time, the loop evolves into a series of 
individual threads and subsequently fades away (Figs. 01 
[5j and Movie 1). 

3.3. Multithermality of the Loop as a function of k 

Since the loop is observed in all six AIA coronal EUV 
bands, it is expected to be mu lti-thermal. We used the 
regul arized inversion method of lHannah fe Kontarl (201 2. 
l2Q13h to obtain the DEMs as a function of n. The results 
are shown in Figs. [6] and [71 

The DEM k reconstruction is performed for each pixel 
of the 1-minute averaged AIA data observed at 13:19 
UT (Fig. 01 second row). Before the DEM reconstruc¬ 
tion, we remo ve the stray light from the AIA data using 
the method of lPoduval et all (j2Q13h . In the DEM recon¬ 
struction, a temperature interval of log(T/K) = 5.5-7.4 
(Alog(T/K) = 0.1) is used. The value of n is assumed 
to be constant throughout the field of view. We note 
that such assumption may not be justified, as individual 
emitting and indeed overlapping structures could in prin¬ 
ciple have different values of n. However, in the absence 
of simultaneous stereoscopic and spectroscopic observa¬ 
tions throughout the field of view, it is not possible to 
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Figure 6. DEM diagnostics from AIA data using the regularized inversion method. The distribution assumed is indicated at each panel. 
Black and green lines stand for the average DEMs within the black and green boxes in Fig.0] Red stands for their difference. Uncertainties 
of the DEMs in each temperature bin are given by the respective horizontal and vertical error bars. 

A color version of this image is available in the online journal. 


obtain the value of ft for each emitting structure in each 
pixel at a given time; an assumption is therefore neces¬ 
sary. By adopting a single value of ft for the entire field 
of view, we demonstrate the feasibility of the regularized 
inversion method to recover the non-Maxwellian DEMs 
from AIA observations using the correspondin g responses 
for the /^-distributions (|Dzifcakova et aO2015l h which (to 
our knowledge) has not been done before. 

The background-subtracted DEMs averaged over the 
small box along the loop (Fig. 0]) are shown in red 
in Fig. [6l DEMs corresponding to the black box are 
shown in black, while the background DEMs correspond¬ 
ing to the green box in Fig. 0]are shown in green. The 
background-subtracted DEMs peak at log(T/K) = 6.2 for 
the Maxwellian distribution and ft = 5. The peak of 
the DEM is shifted to higher T for ft = 2-3, being at 
log(T/K) = 6.3 and 6.4 for ft = 3 and 2, respectively. This 
behaviour of the DEMs is mainly the result of the shifts 
in the ionization equilib rium to higher T, see Fig. [2] and 
iMackoviak et al.l ( 2014 ). The DEMs are multi-thermal, 
with significant amount of emission originating at tem¬ 
peratures lower than the peak of the DEMs, down to 
about log(T/K) = 5.9-6.1. 

The pixel-by-pixel reconstructed EM^(T) are shown 
in Fig. [7] for the five temperature bins where the 
background-subtracted DEM^(T) is the highest (Fig. 
[6|) Note that the width of the temperature bin is 
Alog(T/K) = 0.1. We see that the loop is present in all 
of these temperature bins. At lower temperatures, it is 
relatively well-defined, with the EM^(T) becoming more 
fuzzy with increasing T for all ft. This behaviour comes 
from the emission morphology in the progressively hotter 
AIA bands. 

The DEMs for the Maxwellian distribution and 
ft = 5 contain a spurious high-temperature peak at 


log(T/K) ^ 7.0. This peak is about a factor of ^30 
weaker than the main one and was rep o rted fo r on-disk 
AIA observations also by iDudfk et al.l (j2014cL Sect. 3 
therein). This peak is present also in the background 
DEMs, which may suggest that it is an artifact of the 
method. The peak gets progressively suppressed with 
decreasing ft, until it is absent from the background- 
subtracted DEMs for ft = 2. 

We note that in principle, a combination of AIA fil¬ 
ter ratios permits diagnostics of ft, b ut only for isother¬ 
mal or near-isotherm al structures (jDudfk et al.l 120091: 
iDzifcakova et al.ll2012f h This is because such color-color 
diagram are not monotonic, but contain complicated 
curves. Since the DEMs obtained are significantly multi- 
thermal for all ft studied, it is not possible to use the 
combinations of AIA filter ratios to diagnose ft. Instead, 
the diagnostic of ft has to be performed using ratios of 
individual spectral lines that produce monotonic ratio- 
ratio diagrams (see Sect. 4.3). 

4. HINODE/EIS OBSERVATIONS DURING HOP 226 

The Extrem e -Ultra violet Imaging Spectrometer (EIS, 
iCulhane et al.l [20 071 onboard the Hinode mission 
(jKosugi et al.l 120071 ) provides EUV spectra of the Sun 
in the wavelength ranges 171-212A and 245-29lA with 
a spectral resolution of about 22 mA and spatial resolu¬ 
tion down to 1-2", corresponding to the width of the slit 
chosen. 

On 2013 March 30, EIS was observing portion of AR 
11704 as a part of the Hinode Operation Plan (HOP) 226. 
HOP 226 was originally designed for diagnostics of non- 
Maxwellian ft-distri butions using weak O IV-0 V and 
S X-S XI lines (see IMackoviak et al.l l2013h . Since these 
lines are weak, EIS was rastering an area within the AR 
core using 2" slit and long exposures (60s and 600s). 
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Figure 7. Emission measure EM«(T) at different temperatures T as a function of k. The values of k and log(T/K) are indicated at each 
frame. 

A color version of this image is available in the online journal. 
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Table 1 

Background-subtracted EIS line intensities averaged over the pixels listed. See text for details of the background subtraction. The 
standard deviations cr 10 % and cr 2 0 % are obtained by considering a 10% and 20% calibration uncertainty, respectively. Selfblends are 

indicated, with very weak transitions in parentheses. 


Ion 

A [A] 

selfblending transitions [A] 

Loop (288:317) 

1 010% CO 020%C0 

Loop (300:309) 

1 010% CO 020% (0 

Fe XI 

182.167 

- 

795 

99 

169 

934 

117 

199 

Fe XI 

188.216 

- 

1638 

172 

332 

1947 

204 

394 

Fe XI 

257.554 

257.538, 257.547, 257.558 

398 

45 

82 

414 

47 

86 

Fe XI 

257.772 

257.725 

178 

23 

38 

234 

30 

50 

Fe XII 

186.887 

186.854, (186.931) 

1406 

145 

283 

1498 

154 

302 

Fe XII 

195.119 

195.179, (195.078), (195.221) 

2256 

228 

453 

2506 

254 

503 

Fe XIII 

196.525 

- 

261 

27 

53 

223 

23 

45 

Fe XIII 

202.044 

- 

1346 

153 

279 

1779 

202 

368 

Fe XIII 

203.826 

203.772, 203.795, 203.835 

2591 

270 

524 

2532 

264 

512 


Each 60s and 600s raster contained 10 exposures of the 
2" slit at contiguous positions in Solar X, covering 512" in 
Solar Y (heliocentric co-ordinates). For context, as well 


as density diagnostics, the raster contained additional 
lines, notably several strong Fe XI-Fe XIII lines. 

Upon examination of the data, we find that the 600 s 
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Table 2 

Electron densities diagnosed using individual ratios. 


Line ratio 

h/h 

Loop (288:317) 

log(n e )Maxw. log(n e )^=2 

h/h 

Loop (300:309) 

log(n e )Maxw. log(n e )^=2 

FeXI 182.167/188.215 

Fe XII (186.854+186.887) / (195.119+195.179) 
Fe XIII 196.525 / 202.044 

Fe XIII (203.826 sbl)/202.044 

0.49 

0.62 

0.19 

1.92 

9 45+ u - Ub 

y - yt) -0.08 
y.ou_o ii 

n oc+0.10 

y.oo 0 io 

9 34 + u - 12 
y -°^- 0.20 
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y.55_o 09 

9 20 +0 ‘ 15 
y - ZU -0.15 

9 27 +0 - 13 
-0.14 

0.48 

0.60 

0.13 

1.42 

9 43+ u - Ub 

9 79+ 0 - 10 
‘ y —0.08 

9 12 +0 ' 12 
y -- LZ — 0.12 

y - i 0 -o.io 

9 32 +u ' 12 

y.oz_0 19 

9 76+ 0 - 09 
‘ u - 0.02 

9-02±g;l| 

9.08i°;« 

adopted value 
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)og 10 (AIA 193A [DN s' 1 px' 1 ]) log 10 (l E is l er 9 cm" 2 s "' sr "']) 



Figure 8. A portion of the EIS raster selected for analysis, cor¬ 
responding to the field of view shown by white box in Fig. 0] 
Numbers give the pixel coordinates in the raster. Left: AIA 
193A pseudo-raster corresponding to the EIS FOV. Right : EIS 
intensities of the Fe XII 195.12A+195.18A selfblend. Pixel size is 
2" x T. 

A color version of this image is available in the online journal. 

rasters suffer badly from an accumulation of cosmic rays, 
and that the weak O and S lines do not have sufficient 
intensities for diagnostics. However, the last 60s raster, 
starting at 13:11 UT and progressively rastering in the 
West-East direction, captures a portion of the coronal 
loop in its ninth exposure at 13:19 UT (position 9 in 
Fig. [8]), with the portion of the loop lying directly along 
the EIS slit. Some of the O and S lines proposed for di¬ 
agnostics are again too weak in this 60s exposure. Nev¬ 
ertheless, the strong Fe XI-Fe XIII lines (Table [Tj) are 
sufficient to perform diagnostics of the /^-distributions. 
It is these data that are analyzed here. 

4.1. EIS Data Processing and Calibration 

We first performed coalignment of the raster with the 
AIA data. To do this, we use the Fe XII 195.119A self¬ 
blend (hereafter, “195.119A sbl”) and the AIA 193A data 
uncorrected for s olar differential rotat ion, but with re¬ 
moved stray light (jPoduval et al .1120131 ) . These AIA data 
are used to build an AIA 193A pseudo-raster image by 
extracting the EIS field of view during each exposure 
within the raster, averaging individial AIA 193A frames 
within the duration of each raster exposure, and convolv¬ 


ing with the EIS point-spread function, which is assumed 
to be a G aussian with FWHM of 2" in both X and Y 
directions (jDel Zanna et al.ll2011bl Appendix A therein). 
By comparing this AIA 193A pseudo-raster image with 
the EIS 195.119A sbl, we found shifts in the EIS po¬ 
sitioning of AX = 13" and AT = —9.5" with respect to 
AIA. Furthermore, the EIS slit is found to be rotated 
by +1° with respect to the Solar T. Note that since 
the EIS raster has only 10 exposures in Solar X, the 
EIS raster image cannot be rotated. Rather, the AIA 
data have to be rotated by —1° before constructing the 
AIA 193A pseudo-raster image. Because of the relative 
rotation between the two instruments, as well as for sim¬ 
plicity, the positions of the individual EIS pixels will be 
given in pixel units x and y rather than Solar X and T. 

The EIS data were calibrated using the standard 
eis-prep.pro routine available within the SolarSoft plat¬ 
form running under the Interactive Data Language 
(IDL). During the calib ration, the latest in-flight ra¬ 
diometric calibration of iDel Zannal ( 2013a ) was used. 
The calibrated data contain intensities in physical units 
[erg cm -2 s -1 sr -1 A -1 ] and the 1 -a errors on the ob¬ 
served intensities at each spatial and wavelength pixel. 
The intensities and uncertainties are stored in the sepa¬ 
rate files, with the error file also containing flags for miss¬ 
ing pixels (dusty, hot, warm pixels and pixels affected by 
cosmic rays). We note that usually, these missing pixels 
are interpolated during the calibration procedure. How¬ 
ever, the number of missing pixels in the raster is rela¬ 
tively high, up to 30%. This brings up concerns about 
uncertainties in the fitting of the obs erved intensities at 
each single spatial pixel (|Youndf201Q[ ). 

In this work, the pixels flagged as missing were inter¬ 
polated only for purposes of visualising the raster field of 
view containing the loop (Fig.[8j right). The missing pix¬ 
els were excluded from any further analysis, since they 
cannot always be interpolated with accuracy due to their 
high number. Instead, we rely on intensities obtained by 
averaging along a selected loop segment observed in po¬ 
sition 9 of the EIS slit. The intensity averaging is done 
over y at each wavelength pixel using exclusively pixels 
not flagged as missing. The same is done at position 5 of 
the EIS slit, which we chose to represent the background. 
We selected two loop segments over which the averaging 
is performed. The first one consists of pixels y = 288 to 
317, denoted simply as “loop (288:317)”, and represent¬ 
ing the average loop spectrum. The second one is much 
shorter, consisting of pixels y = 300 to 309, denoted as 
“loop (300:309)”, and representing a much shorter seg¬ 
ment of the loop. We chose to average over longer rather 
than shorter loop segments in order to minimize errors 
from photon noise. We also note that the effective EIS 
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Figure 9. Diagnostics of density using Fe XI, Fe XII, and Fe XIII ratios. Black col or corresponds to th e Maxwellian distribution; red 
stands for k = 2. Different linestyles correspond to different T for each distribution fsee lDudfk et al.l (I2014blj , Figs. 5-7 therein): Full lines 
stand for the temperature corresponding to the peak of the ion abundance, dashed and dot-dashed lines for the temperatures at which 
the relative ion abundance is 1% of its maximum. Observed value of the ratio is plotted as the full cyan line. Uncertainties including 
calibration uncertainties are plotted for illustration. 

A color version of this image is available in the online journal. 
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Fe XI 188.216A / 

/ Fe XI (257.538A + 257.547A 4- 257.554A + 257.558A) 



Fe XI 188.216A / 

/ Fe XI (257.538A + 257.547A 4 257.554A ■ 


257.558A) 



Figure 10. Diagnostics of k, involving the Fe XI 257.554A selfblend. Different colors correspond to different k, linestyles denote electron 
density. Diamonds and asterisks denote the predicted ratios based on DEM analysis. 

A color version of this image is available in the online journal. 

resolution in y is only about 3 // -4 // because of the instru- ment point-spread function. 









































































13 








Figure 11. Diagnostics of k involving the Fe XI 257.772A selfblend. Different colors correspond to different k, linestyles denote electron 
density. Diamonds and asterisks denote the predicted ratios based on DEM analysis. 

A color version of this image is available in the online journal. 
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Subsequently, the average spectrum of the background 
at position 5 is subtracted from the average loop spec¬ 
trum. The resulting background-subtracted loop spec¬ 
trum is fitted with Gaussian line profiles and locally lin¬ 
ear continua, taking into account known bl e nds and self- 
blends (e.g., iDel Zanna] I 2010 L 1201 lL 120121: iDndfk et al.l 
l2014bl) . The intensities obtained are listed in Table 
[TJ The uncertainties on these intensities are obtained 
by error propagation of the uncertainties of the photon 
noise and dark-current subtraction that are the output of 
the eis-prep.pro. These uncertainties are then added in 
quadrature with the uncertainty of the EIS radiomet¬ 
ric calibration. We consider two values of the radio- 
metric calibration uncertainty, 10% and 20%. The 20% 
uncertainty is the un certaint y of the ground cal ibration 
(jCnlhane et al.l l2007f ). while IWang et all (j 2011 h argues 
that the in-flight uncertainty is smaller, about 10%. The 
final uncertainties on the line intensities are denoted cr 10 % 
and <7 2 o% 5 respectively, and are also listed in Table [lj 

We note that we employed the revised radiomet¬ 
ric calibration of IDel Zannal (12013af) r ather than the 
ground calibration (iCulhane et al.l 2007 ) or the calibra¬ 
tion of lWarren et alJ ~ (j2014f ). IDel Zannal (j2013al ) showed 
that significant departures from the EIS ground cal¬ 
ibration occurred over time, especially with the lines 
in the long-wavelength channel of EIS being underes¬ 
timated by about a factor of two for observations af¬ 
ter 2010. Line ratios were used to obt ain a calibra¬ 
tion c orrected for decrease of sensitivity. I Warren et al.l 
(|2014f ) obtained similar results, comparing the EIS ra¬ 
diances from the whole Sun with the SDO /EVE irradi- 
ances together with the quiet-Sun DEM modeling. How¬ 
ever, the EVE calibration can be an additional source 
of uncertainty. Furthermore , the q uiet-Sun DEM mod¬ 
elling used by Warrer Tet al.l ( 20141) can in itself be sen¬ 
sitive to ft (jMackoviak et al.ll2014f ). which could entan¬ 
gle the calibra tion to diagnostics and overcomplicate the 
problem. The IDel Zannal (j2013al ) calibration is sufficient 
for our purposes, as the diagnostics of electron density 
(Sect. 14.21 ) as well as the d iagnostics of ft (Sect . 14.31) 
rely only on line rati os (Dzifcakova & K nlinovdl 120101: 
iMackoviak et al .1120131: iDudfk et al.ll2014b[ ). and thus do 
not require absolute calibration. 

We also note that the decrease of sensitivity of the 
long-wavelength channel of EIS by about a factor of 2 
is clearly significant for our diagnostic purposes. This 
is because the diagnostics of ft in Sect. 14.31 involve lines 
from both the short-wavelen gth an d long - wavele ngth EIS 
channels. We note that the IDel Zannal (|2013aD calibra¬ 
tion applies only to observations up to September 2012 , 
therefore we have assumed that no further degradation 
occurred since. If the long-wavelength channel sensitiv¬ 
ity further degraded significantly, the radiances of the 
long-wavelength channel would be underestimated, and 
the ratios shown in Figs. [T01TL21 would decrease. 

4.2. Density Diagnostics 

Diagnostics of electron density are a neccessary pre¬ 
requisite to diagnostics of ft (Sect. 14.3|) . The Fe 
lines observed by EIS (Table [Tj) contain several com¬ 
binations of lin e s sen s itive to electron d e nsity (e.g., 
Watanabe et al l2009t lYou ng et al. 120091: Del Zanna] 


2010 , 201 It IDel Zanna et al.ll 2012 bl: IDndfk et al 
The density-sensitive line ratios are listed in Table 


where the diagnosed densities are also listed. 

The diagnostics of density are performed using 
the m ethod developed first by iDzifcakova fe Kulinoval 
(1201 Of). The theore tical line ratios were calculated by 
IDndfk et al.l (|2014bf ) using the latest available atomic 
data (see also Sect. 12 . 21 ) . This density diagnostics 
method employs ratios of lines arising in the same ion. 
The sensitivity of these ratios to n e cannot be signifi¬ 
cantly influenced by their sensitivity to T and ft, other¬ 
wise the density diagnostics would be precluded. Typi¬ 
cally, it was found that density-sensitive ratios commonly 
used for the Maxwellian distribution can be used under 
the assumption of ft-distributions in the same manner, al¬ 
beit the resulting densities can be up to 0.1 dex lower for 
ft = 2 (red color in Fig. (9j This uncertainty due to the un¬ 
known value of ft is however comparable or smaller than 
the uncer tainty due to the dependence of individual ra¬ 
tios on T (Dzifcakova & Kulinova 2010; M ackoviak et al.l 
l2013t IDndfk et al.ll2014a l lbh . The overall uncertainty due 
to the a-priori unknown values of T and ft can be in 
some instances as large as 0.4 dex in log(n e /cm -3 ) (e.g., 
Fig. |9j). 

In principle, calibration uncertainties compound the 
uncertainty of the diagnosed density. However, the EIS 
density-sensitive lines belonging to the same ion have 
typically similar wavelengths, and are all observed within 
the same EIS channel. It is unlikely for the correspond¬ 
ing line ratios to be severely affected by calibration un¬ 
certainties. Therefore, we chose to ignore the calibra¬ 
tion uncertainties when diagnosing the electron density. 
Nevertheless, for illustration, the values of individual 
ratios deriv ed fr om the background-subtracted observa¬ 
tions (Sect. 14.ip are plotted in Fig. [9] together with their 
uncertainties. 

The densities diagnosed using the Fe XI and Fe XIII ra¬ 
tios are consistent within their respective errors (Fig. [9j 
Table [2]). However, the densities diagnosed using Fe XII 
are higher by about 0.5 dex compared to those diag¬ 
nosed using Fe XIII. I.e., given the densities diagnosed 
using Fe XI and Fe XIII, the calibrated intensities of 
the Fe XII 186.887A selfblend seem to be too high com¬ 
pared to the intensities of the Fe XII 195.119A sbl by 
about 2 0-30%. The reason for t his discrepancy is un¬ 
known. IDel Zanna et al.l ( 2012 b ) showed that the new 
atomic model for Fe XII provides significantly lower 
densities than the previous model, by about 0.4 dex. 
The previous model consistently provided values much 
high er than those obtained from Fe XI and Fe XIII (see, 
e.g.. lYoung et al.ll2009[ ). The densities obtained for the 
background-subtracted intensities from Fe X II are much 
reduced using the new atomic data of IDel Zanna et ahl 
(l 2012 bl) . but still higher by about 0.5 dex than those ob¬ 
tained here from Fe XI and Fe XIII. The discrepancy is 
lower when background-subtraction is ignored. The rea¬ 
son is unclear despite the fact t hat we use here the same 
ato mic data an d calibra tion as IDel Zanna et al.l (| 2012 bf ) 
and IDel Zannal (|2013af ). Nevertheless, the intensity of 
the 186.887A line is somewhat susp ect since the shape 
of the EIS effective area curve in the IDel Zannal (j2013al ) 
calibration differs from that of the ground calibration 
(see Fig. 8 top therein), or from that of I Warren et al.l 
(| 201 4l . Fig. 7 therein). The differences are of the or¬ 
der of several tens of per cent, which is comparable to 
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Figure 12. Diagnostics of k involving the Fe XI 257.772A selfblend and Fe XII / Fe XI ratios using lines only from the EIS short-wavelength 
channel. Different colors correspond to different k, linestyles denote electron density. Diamonds and asterisks denote the predicted ratios 
based on DEM analysis. 

A color version of this image is available in the online journal. 
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the 20% calibration uncertainty. We note that chang¬ 
ing the Fe XII 186.887Asbl /195.119Asbl ratio by about 
20% would bring the diagnosed densities into consistency 
with the Fe XI and Fe XIII ones (Fig. [9]). Furthermore, if 
the 20% calibration uncertainties are taken into account, 
the observed Fe XII ratio still yields densities consistent 
with the other ratios, since the lower violet line in Fig. 
[9] yields densities of log(n e /cm -3 ) = 9.2-9.4. 

Additionally, because of the strong overlap of the 
Fe XI-Fe XIII relative ion abundances (Fig. [2j) at tem¬ 
peratures corresponding to the peaks of the background- 
subtracted DEMs for all ft (Fig. [6]), it is unlikely that 
the diagnosed loop is a multi-density one. If it were, the 
Fe XIII density-sensitive ratios, which have have stronger 
density sensitivity at log(n e /cm -3 ) = 9.5, should be bi¬ 
ased towards higher densities more than the Fe XII ratio, 
whose increase with n e is weaker. 

For these reasons, we adopt the densities diagnosed 
using Fe XI and Fe XIII for further analysis. For 
the average loop spectrum, the diagnosed value is 
log(n e /cm -3 ) =9.0-9.5, while for the loop spectrum av¬ 
eraged through pixels y = 300-309 we obtain the value of 
log(n e /cm -3 ) = 8.8-9.6. We note that these values are 
conservative, as they are chosen to contain all the den¬ 
sities diagnosed using Fe XI and Fe XIII together with 
their respective errors (Table EJ). 


4.3. Diagnostics of T and ft 
4.3.1. Ratio-Ratio Diagrams 


Having obtained constraints on the electron density, we 
can next diagnose ft. A principal constraint is that the 
diagnostics of ft have to be performed using lines that are 
also sensitive to temperature. This comes from the na¬ 
ture of the task, as both the excitation as well as ioniza¬ 
tion and recombination rates are a function of both T and 
ft. Typically, line ratios sensitive to ft involv e lines pop¬ 
ulated by different parts of the distribution (|Dzifcakoval 
l2006at h i.e., either lines with different behaviour of the 
excitation cross-section with E; or lines with different 
excitation thresholds A Eji, i.e., lines formed at different 
wavelengths; or both. Such line ratios will always be sen¬ 
sitive to temperature as well. Line ratios sensitive to ft 
are then be combined with the line ratios involving lines 
from neighbouring ions, that are strongly sensitive to T, 
but also sensitive to ft. 

Furthermore, if the relative level population de¬ 
pends on density, /^-sensitive line ratios belonging to 
the same ion typically have smaller sensi tivity to ft 
than t o n e , although except i ons occur dDzifcakova 
2006a; [D zifcakova & Kulinoval 120101 : IMackoviak et al. 


20131 : Dndfk et al. 2014b). The diagnosed electron den¬ 


sity is then used to co ntrain the diagnostics using 
the ratio-ratio diagra ms (Dzifcak ova fe Kulinoval [20101: 
IMackoviak et al1l2013[ ). 

This is the approach we chose here. From the 
lines listed in Table [lj Fe XI is the only ion offering 
strongly temperature-sensitive line ratios. There are 
four temperature-sensitive ratios of Fe XI lines, each in¬ 
volving a line from the EIS short-wavelength channel 


(182.167A or 188.216A) and a line from the EIS long- 
wavelength channel (257.554A sbl or 257.772A sbl, see 
Table [Tj) . These Fe XI line ratios are combined with 
ratios involving a Fe XII line. Both the selfblends at 


186.887A and 195.119A are used. This is because of the 
possib le ca libration problems involving the 186.887A sbl 
(Sect. 14.2p . The various combinations are shown in Figs. 
M ID and [121 The observed line ratios together with 
their uncertainties are shown as large crosses where the 
azure and violet errorbars are calculated using the cr 10 % 
and a 2 o% uncertainties of the individual lines involved, 
respectively (see Table [Tj) . The calibration uncertainty 
is included since various ratios involve lines observed in 
both EIS channels. 

The ratio-ratio diagrams involving the 257.554A sbl 
(Fig. ITU|) show strong sensitivity to T and sensitivity to 
ft that is of the order of ^20%. However, the large cal¬ 
ibration uncertainties together with the uncertainty in 
the diagnosed log(n e /cm -3 ) prevent positive diagnostics 
using some of these diagrams for both loop segments. 
This is because the l-cr 2 o% errorbar intersects all curves 
for different ft for at least one of the diagnosed limits 
on log(n e /cm -3 ). The situation is worse for the loop 
(300:309), where the observed ratios are closer and in¬ 
tersect more curves for different ft on the ratio-ratio dia¬ 
grams (Fig. [lOj left ) than in the case of the average loop 
spectrum (loop (288:317), Fig. [TUI right). For the aver¬ 
age loop spectrum, we can diagnose ft < 5 using the ratios 
involving Fe XI 188.216A / (Fe XI 257.554A sbl). This is 
because for all ratios in Fig. [TO] right , the red and orange 
ratio-ratio curves corresponding to ft = 2 and 3 intersect 
all the violet errorbars, while others do not. Neverthe¬ 
less, the green ratio-ratio curve (ft = 5) is in close vicinity 
of the errorbar in the top right panel. 

The situation is much better when using ratio-ratio di¬ 
agrams involving the Fe XI 257.772A selfblend (Fig. HI]) . 
Compared to the spread of the curves for individual ft and 
log(n e /cm -3 ), the uncertainty of the observed ratios is 
smaller, and the observed ratios are located further away 
from the curves. This permits diagnostics of ft <3, a 
strongly non-Maxwellian distribution. We note that this 
result is the same whether the Fe XII 186.887A selfblend 
or the 195.119A selfblend is used. 

To confirm this diagnostic, in Fig. [12] we substitute 
the Fe XI 257.772 A sbl by other Fe XI lines in the 
Fe XII / Fe XI ratio. Again, ft < 3 is found independently 
of the combination of lines used (Fig. [12]), with majority 
of the ratio-ratio diagrams indicating ft <2. 

We note that the results of the diagnostics of ft are not 
dependent on which of the two Fe XII lines is used. Both 
the Fe XII 186.887A and 195.119A selfblends yield similar 
results in terms of ft (Figs. [T01TL2]) . This allows us to 
establish confidence in using the Fe XII lines to diagnose 
ft despite the possible calibration problems mentioned in 
Sect. 14.21 

Finally, we point o ut th at the loop evolution as ob¬ 
served by AIA (Sect. I3.2j) is unlikely to strongly affect 
the diagnostics of ft. This is b ecause the AIA 1 93 A in¬ 
tensities (dominated by Fe XII; iDel Zannall2013bh do not 
change by more than 1-2% during the EIS observations 
at 13:19UT (Fig. 0 top). 

4.3.2. Influence of the DEM on Diagnostics of k 

Originally, the ratio-ratio diagrams were developed 
for simultaneous diagnostics of T and ft based on the 
assumption that the observed structure is isotherma l 
([Dzifchkova fe Kulinoval 120101 : IMackoviak et al.l I2013D . 
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Figure 13. Same as Fig. 1111 but with atomic data corresponding to CHIANTI 7.1 and the KAPPA package. 
A color version of this image is available in the online journal. 
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The results presented for the transient loop in Sect. 13.31 
however suggest that the loop is multi-thermal indepen¬ 
dently of the value of ft. Therefore, we investigated the 
influence of DEM on the diagnostics of ft. 

Unfortunately, the EIS raster obtained during HOP 
226 does not contain enough strong lines for determi¬ 
nation and constraining the DEM, especially through¬ 
out the entire log(T/K) = 5.5 - 7.0 range. Therefore, 
we use the AIA d ata to per f orm th e DEM diagnos¬ 
tics and note that iDel Zannal (2013b) showed that the 
AIA observations can be used to predict the EIS ra¬ 
diances to within the calibration uncertainties of both 
instruments. Although these authors used preferen¬ 
tially a DEM recon s truct ion technique different from the 
Han nah fe Kontarl (120121) one that is used here (Sect. 
E3D, the lHannah fe Kontarl (j2012h technique was also 
tested by Del Zanna (|2013bh and a reasonable agreement 
was found. 

First, we produce an AIA pseudo-raster for each AIA 
band similar to that shown in Fig. [8j left. Subse¬ 
quently, the DEM is derived using the technique of 
lHannah fe Kontarl (|2012l . l2013h similarly as in Sect. 13.31 
for each pixel of these AIA pseudo-rasters. We perform 
the averaging of the DEMs the same way as for the EIS 
raster (i.e., over pixels at position 9 corresponding to the 
loop, y = 288-317 and 300-309), subtract the background 
at position 5, and then use these background-subtracted 
DEMs to predict the intensities of individual EIS Fe XI- 
Fe XIII lines. 

The predicted ratios of individual lines are shown as 
a function of ft on each of the ratio-ratio diagrams in 
Figs. [TQTH21 Diamonds are used for the diagnosed 
lower limits on log(n e /cm -3 ) (Sect. 14.2|h asterisks for 
the upper limits. Generally, the predicted ratios for 
each ft are located close to the corresponding curves, 
indicating that these curves can still be used to indi¬ 
cate the value of ft in the observed plasma even if this 
plasma is multi-thermal. Furthermore, we find that with 
decreasing ft, the predicted ratios converge on the ob¬ 
served values in all cases, indicating that the plasma 
is strongly non-Maxwellian with ft <2. This is so even 
for the cases when the isothermal ratio-ratio diagrams 
cannot be used to constrain the value of ft, such as 
the Fe XI 188.216A / (Fe XI 257.554A sbl) - (Fe XII 
186.887A sbl) / (Fe XI 257.554A sbl) in Fig. [lOj middle 
row , or the similar combination with Fe XII 195.119A sbl 
in Fig. ca bottom row. 

5. ATOMIC DATA UNCERTAINTIES 

The atomic datasets for astrophysical spectroscopy are 
always incomplete, since they contain only a finite num¬ 
ber of energy levels and the corresponding transitions. 
Therefore, we investigated the influence of the atomic 
data uncertainties on the analysis presented here. To do 
that, we repeated the analysis presented in Sect. [4] us¬ 
ing older at omic data that are available in the C HIANTI, 
version 7.1 (jDere et al.lll997l:lLandi et al.ll2013[ ). The cal¬ 
culations for the Maxwellian distribution were performed 
using CHIANTI v7.1 and the corresponding calculations 
for the ft -distributions were perfo rmed using the KAPPA 
package (Dzifcak ova et al1l2015[ ). 

The density diagnostics performed in Sect. 14.21 using 
Fe XI and Fe XIII ratios remain valid, since these density- 


sensitive ratios do not chang e appreciably (see Figs. 5 
and 7 in iDudfk et al.l l2014b l even if the individual line 
intensities change by up to ~20%. Therefore, the con¬ 
servative densities adopted (Table |2j) are kept. 

Examples of the diagnostics of ft, including the ef¬ 
fect of DEM using the older atomic data, are presented 
in Fig. [13j which shows ratios involving the Fe XI 
257.772A selfblend. These diagrams correspond to those 
shown in Fig. [IT] for the newest atomic data. From 
Fig. d3] we see that the diagnostics of ft < 2 remain valid 
even if older atomic data are used. 

We note that the newest atomic data of IDel Zannal 
(I2Q1QT) . IDel Zannal (l2011lh Del Zanna et al.l (j2012bf h and 
IDel Zanna fe Storevl (|2Q13h used in Sects. [2H4] represent 
significant improvement over the previous ones, espe¬ 
cially in the case of Fe XII, where differe nces of up to 60% 
in intensities of key lines were found (jDel Zanna et al.l 
l2012bf h including the 195.119A selfblend. We find sim¬ 
ilar increases found for all ft. The Fe XI line intensi¬ 
ties are found to be increased for all ft by about 10- 
20%, although details depend on ft and the particular 
line. These increases in line intensities are due to the 
increased contribution from resonances as we l l as ca s¬ 
cading from the n = 4 lev el (De l Zanna et al.l l2012allbl : 
IDel Zanna fe Storevl 120131 ) and cause the change in the 
theoretical ratio-ratio diagrams (compare Fig. [lT] with 
Fig. fl3l) . 

Finally, we note that the newest atomic data used 
in Sects. EHU do not include contributions from the 
n ^ 5 levels. For Fe XI, the contributions of cascading 
from ft > 5 to Fe XI line intensities are about 10-20% 
(|Del Zanna fe Storevll2013 ). For the transitions from the 
3s 2 3p 2 3d configuration in Fe XII (i.e., the 186.887A and 
195.119A lines), the missing contributi ons from casca d¬ 
ing is smaller, of the order of 5% only (jDel Zanna et al.l 
l2012bh . Including these contributions would move the di¬ 
agnostic diagrams (Fig. [T0l4T2l) down in the ^/-direction. 
Therefore, more complete atomic data are unlikely to 
change the results of the diagnostics of ft. 

6. SUMMARY 

We reported on imaging and spectroscopic observa¬ 
tions of a transient coronal loop observed within the core 
of AR 11704 on 2013 March 30. The loop reappeared in 
the same location as the already faded flare arcade of the 
B8.9-class microflare peaking about three hours earlier. 
The transient loop persisted for nearly two hours, during 
which it evolved into a series of individual threads. By 
examining the AIA data, we found no signatures of hot, 
flare-like emission being associated with the loop. These 
results were confirmed using the DE M reconstruction b y 
the re gular ized inversion met h od of lHannah fe Kontarl 
(|2012f ) and lHannah fe Kontarl (|2013f h showing the rela¬ 
tive absence of hot plasma except a spurious peak that 
is also present in the background. This method was em¬ 
ployed in conjunction with the AIA r esponses emission 
calcul ated for the ft-distributions by iDzifcakova et al.l 
(|2015f ). We found that the loop is multi-thermal for all 
ft considered, with the DEMs peaking at log(T/K) = 6.3 
for the Maxwellian distribution, and at 6.5 for ft = 2. The 
spurious high-temperature peak becomes less prominent 
with decreasing ft, and it disappears for ft = 2. 

We analyzed the spectroscopic data obtained by Hin- 
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ode/ElS in order to perform diagnostics of density and 
ft. Several density-sensitive ratios of Fe XI, Fe XII and 
Fe XIII were used in conjunction with the latest EIS 
calibration and the latest available atomic data. We 
find consistency between the densities diagnosed from 
the background-subtracted Fe XI and Fe XIII intensi¬ 
ties, with the Fe XII intensities yielding higher densities. 
However, the density-diagnostics using Fe XII may still 
be consistent with the Fe XI and Fe XIII ones if the 
EIS calibration uncertainty is taken into account. Using 
Fe XI and Fe XIII yields a conservative estimate of elec¬ 
tron density of the order of log(n e /cm -3 ) = 8.8-9.6 for 
a loop segment 10" long, and log(n e /cm -3 ) =9.0-9.5 for 
the average spectrum of the loop along position 9 of the 
EIS slit. 

This diagnosed electron density is then used to con¬ 
strain the diagnostic of ft. To do that, we use the 
temperature-sensitive ratios involving Fe XI lines from 
both EIS detectors in combination with the ratios involv¬ 
ing Fe XII and Fe XI lines. We found that ratios involving 
the Fe XI 257.554A selfblend can preclude a succesful di¬ 
agnostic due to the calibration uncertainty being larger 
than the spread of the curves for individual ft. However, 
all ratios involving the Fe XI 257.772A selfblend together 
with other Fe XI and Fe XII lines consistently yield ft < 2, 
i.e., an extremely non-Maxwellian situation. 

We next studied the influence of the plasma multither- 
mality on the diagnostics of ft. Since the EIS data are 
insufficient to obtain the DEMs, we obtained the DEMs 
from the AIA data. These DEMs were used to predict 
the EIS line intensities as a function of ft. We found that, 
with decreasing ft, all ratios of the predicted Fe line in¬ 
tensities converge on the observed values. These results 
confirm the diagnosed value of ft = 2 and provide first 
quantitative description of the non-Maxwellian distribu¬ 
tion of electron energies in a coronal loop. 

We note that the transient loop studied here is not a 
typical coronal loop. This is due to its rather low tem¬ 
perature compared to typical active region cores, the on¬ 
going magnetic reconnection as evidenced by the bright- 
enings and jetting activity near its right footpoint, as 
well as its appearance in the same spatial location as the 
previous B8.9-class microflare. Nevertheless, the results 
presented here demonstrate the viability of the methods 
for diagnostics of ft. Such methods should be applied in 
the future on the analysis of typical, well-defined coro¬ 
nal loops to search for possible signatures of impulsive 
heating. 

Finally, we note that the 20% calibration uncertainty, 
typical of the EUV spectroscopic instrumentation, can 
represent a severe limitation to the diagnostics of non- 
Maxwellian distribution. Such large calibration ucer- 
tainty contributes to the possible mis-interpretation of 
the observations, as it is always possible to obtain some 
temperatures and DEMs, if a diagnostic of ft is not per¬ 
formed. Decreasing the calibration uncertainty to about 
10 % would significantly contribute to enabling the diag¬ 
nostics of non-Maxwellian distributions from EUV spec¬ 
tra. 
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